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ABSTRACT 

The Hubble parameter H^ is still not very well measured. Although the Hubble 
Key Project, Chandra and WMAP gave good estimates, the uncertainties remain 
quite large. In this brief report, we suggest an original and independent method 
to derive a lower limit on Hq using the absorption of very high energy gamma- 
rays by the cosmic infrared background. With conservative a hypothesis, we obtain 
Hq > 74 km.s^^.Mpc"^ at the 68% confidence level, which favors the upper end of 
the intervals allowed by dedicated experiments. 



Key words: cosmological parameters - diffuse radiation. 



INTRODUCTION 

Gamma-rays in the 10^^ — 10^^ eV range are observed from 
distant Active Galactic Nuclei (AGN). Those very high- 
energy photons will be absorbed by the Cosmic Infrared 
Background (GIB) due to electron-positron pair production. 
In the higher energy range (typically around 10-20 TeV), 
the absorption is expected to be so high that the difference 
between the source spectrum and the measured spectrum 
becomes very large. As the GIB spectral distribution is 
now quite correctly known (at least with lower limits) it 
is possible to compute the intrinsic spectrum of a given 
' AGN as a function of the integrated GIB density crossed 
along the line-of-sight through an unfolding procedure. The 
absorption being proportional to e^'^", it is then possible 
to derive a sensitive lower limit on the Hubble parameter 
so as to exclude unfolded spectra which would be explicitly 
unphysical. The idea is very simple : the smaller the 
Hubble constant, the larger the distance to the source, the 
larger the number of GIB photons crossed, and the higher 
the unfolding coefficients. As the unfolded {i.e. intrinsic) 
spectrum shape will inevitably become unacceptable for low 
enough Hubble parameters, this approach allows to derive 
a lower limit on Hq. The obtained bound is relevant when 
compared with current estimates, although the main aim of 
this brief article is to open this field and outline the scheme 
of the method. In the first section, the latest measurements 
of the GIB density are reviewed, including some recently 
revised analysis. In the second section, the gamma-ray data 
relevant for this study are presented together with the 
associated uncertainties. In the third section, the physical 
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processes involved in the gamma-ray absorption mechanism 
are explained and the basic features of the proposed method 
are underlined. The fourth section deals with the details 
of the unfolding procedure and with the physical criteria 
used to reject "unphysical" unfolded spectra. In the fifth 
section, the Monte-Carlo method used to derive statistically 
meaningful bounds is described and the results are given. 
Finally, several developments that could be expected in the 
future arc outlined in the last section. 



1 THE COSMIC INFRARED BACKGROUND 
DENSITY 

Our understanding of the early epochs of galaxies has 
recently increased thanks to the observational evidences 
provided by UV/Visible/Near-IR, far-IR and submillimeter 
surveys of high-redshift objects. In a consistent scenario, 
galaxy formation and evolution can also be constrained 
by the background radiation which is produced by the 
line-of-sight accumulation of all extragalactic sources. 
The Cosmic Infrared Background (GIB) is basically the 
relic emission at wavelengths between a few microns and 
millimeters of the formation and evolution of galaxie s of a ll 
type and star-f ormin g systems (see, e. g., Pu get et al. (l99a), 
Hause r et al. (119981 ), Lagache et al. (| 19991 ). G ispert et al. 
(I 2OOO I), Hauser & Dwek |200il) and Kashlinsky (|2005l )). The 
near-IR GIB arises mainly from the stellar component of 
galaxies and probes their evolution at early times. The mid- 
and far-IR GIB originates from dusty galaxies reprocessing 
stellar light and other energetic output. 
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Observationally, the CIB is difficult to distinguish from 
the generally brighter foregrounds contributed by the local 
matter within the solar system, the stars and the interstellar 
medium of the Galaxy. However, the situation has dramat- 
ically improved in the last decade. In this study, we have 
used the most accurate and "up-to-date" esti mates of the 
CIB density, as described, e.g., in Dole et al. l|2006l ') with 
linear interpolations (in log-log scale) between the measure- 
ments. The doubtful measurements have deliberately been 
ignored. The number of points is now large enough to make 
this approach quite accurate. The measurements used in our 
analysis are the following: 

• .3 pm : 10±5 nW.m~ ^.sr~^ , from Bernstein et al. 
l|2002h corrected by Matilla (|2003l '). The estimate is based 
on the Hubble Space Telescope (HST) broadband CCD pho- 
tometry and the 2.5 m du Pont telescope from Lac Cam- 
panas Observatory (LCO). 

• 1.25 urn : 28±15nW.m"^sr~\ from Wright l|200lD and 
Cambresy et al. (|200ll '). The estimate relies on the Diffuse 
Infrared Background Experiment (DIRBE) and 2 Micron 
All Sky Survey (2MASS) data. 

• 2.2 /im : 22.4±6 nW.m'^.sr" 
l|2000l ) using DIRBE and Lick data. 



3.3 /im : 11.1±3.3 nW.m~^sr- 



, from Gorjian et al. 
from the same analy- 



• 15 ^m : 3.3±1.3 nW.m~^sr"\ from Elbaz et al. (|l998l ). 
This lower limit is based on galaxy counts with the Infrared 
Space Observatory Camera (ISOCAM). 

• 2 4 fj,m : 2.7±1.1 nW.m~^.sr~^, from Papovich et al. 
(1200J). The limit comes from galaxy counts with the Multi- 
band Imaging Photometer for Spitzer (MIPS). 

• 70 ^m : 7.1±1.0 nW.m"^sr"\ from Dole et al. (|2006h . 
based on MIPS data with more than 19000 24 /im sources 
stacked with ^24 > 60 fi Jy. 

• 160 fj,ia : 13.4±1.7 nW.m^^.sr^^, from the same anal- 
ysis. 

It is important to underline that most of those esti- 
mates are lower limits. As the constraints derived on the 
Hubble parameter in the next sections would have only 
been weakened by a lower CIB density, this makes the 
following analysis conservative. The Spitzer observatory 
data (Werner et al. (20045)) allowed for deep and wide 
area surveys, in particular at 24, 70 and 160 fim using the 
MIPS (Rieke et al. (|2004l ')') which is a keypoint for this 
analysis. Those point s, used here with the refined a nalysi s 
of Dole et al. (120061 ') . improving over Dole et al. (|2004h . 
are very relevant as they (at least two of them) lie within 
the energy range where the interactions with gamma-rays 
are large. In the following, we therefore consider the CIB 
density as known (within the uncertainties) and use it as 
an "absorbing material" for the gamma-rays. The higher 
the distance to the source (i.e. the lower the Hubble 
parameter), the higher the absorption. As the transmitted 
flux depends exponentially on the depth of the absorber, 
the approach is very sensitive. 

Fig. [1] displays the data points (filled circles) used for 
this analysis together with lower limits (triangles) and es- 
timates not taken into account as being affected by higher 
uncertainties (open circles). The Monte-Carlo method de- 
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Figure 1. Cosmic infrared background energy density as a func- 
tion of the wavelength. The filled circles are the data selected for 
this analysis. Open circles are tentative detections ignored in this 
article and triangles are lower limits. The intervals labeled "4 TeV 
sensitivity" and " 17 TeV sensitivity" correspond to 90% of the in- 
teractions of, respectively, 4 TeV and 17 TeV gamma-rays taking 
into account the cross section convolved with the CIB spectrum. 
The shaded region corresponds to the \a allowed values of the 
CIB density, as included in the Monte-Carlo simulation. 



scribed in the next sections samples CIB realizations within 
the shaded region. 



2 ACTIVE GALACTIC NUCLEUS MRK 501 

It has been suggested long ago that observations of the TeV 
spectrum of extragalactic sources can be a powerful tool 
to constrain the CIB spectrum, especially around 10 /xm, 
where the Extragalactic Background (EB) constraints were 
quite weak. As a matter of fact, TeV gamma rays propagat- 
ing in the intergalactic medium undergo absorption through 
electron-positron pair production on CIB photons. Using 
sources Mrk 421 and Mrk 501, meaningful upper limits 
have been established in different papers. More recently, 
interesting constraints were obtained by Aharonian et al. 
(200q). They are based on the observations performed 
with the very effec tive ne w network of telescopes H.E.S.S. 
(Aharonian et al. l|2005al )). In this article, we reverse the 
argument and use reliable estimates of the CIB density to 
obtain a new constraint on the Hubble parameter. T he ide a 
has already been pointed out by Salamon et al. ([199J). 
However, Mrk 501 data were not available at that time 
and no result could be derived. Furthermore, the method 
suggested by Salamon et al. is far less conservative as it 
assumes that any roUoff in the AGN spectrum should be 
due to the CIB absorption. This is a very strong hypothesis 
as many different effects could (and, even, should) lead to 
a roUoff intrinsic to the source. This could, among other 
processes, be due to self-absorption within the AGN, to 
the Klein-Nishina effect or simply to the limiting energy 
available for Fermi acceleration to take place in a finite 
volume over a finite amount of time. The bound derived in 
this article does not rely on this "CIB rollofF' assumption 
and is therefore much more robust. But it can only lead to 
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a lower limit and not to a measurement. 

The TeV source Mrk 501 is the second closest X- 
ray selected BL Lac object after Mrk 421, with a redshift 
z = 0.034. During the 1997 outburst which lasted several 
months, Mrk 501 w as ob served intensively in X-ray s (Bep - 
poSAX: Plan ei al. (jl998t ). RXTE: Lamer & Wagn er (| 19981 )1 
and TeV 7- rays (Whipple: Catanese et al. (|l997), Samuel- 
son et al. (|l998l '): HEGRA: Ah aronia n et al. (1997), Tele- 
scope Array : Hayashida et al. (|l998l l: CAT: Djannati-Atai' 
et al. (|l999h '). The exceptional 1997 April 16 flare was ob- 
served by BeppoSAX and low-energy threshold (~ 300 GeV) 
Whipple and CAT atmospheric Cherenkov telescopes, al- 
lowing the derivation of the energy spectrum with a good 
accuracy in a broad dynamical range. These unique data 
initiated interesting efforts to set meaningful upper limits 
on the CIB flux ( see e.g. Biller et al. (1998), Stanev & 
Franc eschini ()l998t ). Stecker & De Jager (|l998l ). Aharonian 
et g^. (|l999bl ). Copp i fc A haronian (|l999t ). Konopelko et al. 
j 19991 ). Guy et al. l|2000l ). Renauh et al. (|200ll ')'). More re- 
cently, other AGNs have been detected, especially with the 
H.E.S.S. network (see, e.g., Aharonian et. al. (200513) for a 
review). As far as our study is concerned, a promising op- 
portunity is the measurement of the spectrum of the distant 
blazar lESllOl-232 (Aharonian et al. l|2007l ')'). With a red- 
shift z = 0.186, this is a very good candidate to maximize 
the interaction with the CIB. The spectrum has been mea- 
sured by H.E.S.S. up to ~ 2 - 3 TeV. In fact, the optical 
depth for gamma-rays is roughly the same at this energy 
for z = 0.186 as at ~ 15 - 20 TeV for the lower Mrk 501 
redshift z = 0.034. The analysis hereafter described has also 
been performed for both Mrk 501 and lESllOl-232. How- 
ever, as the results are less stringent for the latter, we will 
mainly focus on Mrk 501. Furthermore, the very wide range 
of measured energies (between ~ 400 GeV and ~ 20 TeV) 
for Mrk 501 allows for a much more reliable analysis. In this 
article, CAT and HEGRA data ar e used from 400 GeV to 
21 TeV (Aharonian et al. (|l999bl ') and Guy et al. (|2000l ')'). 
The higher energy points, i.e. the more relevant ones for 
this analysis, come from a time averaging of the HEGRA 
spectra based on more than 38000 detected 7-rays. This 
very high number of data points combined with the 20% 
energy resolution of the instrument resulted in a very ac- 
curate measurement of the spectrum above 10 TeV which 
remains unequaled. 



3 INTERACTIONS OF GAMMA-RAYS WITH 
THE INFRARED BACKGROUND 

The influence of low energy photons in the Universe on 
the propagation of Very High Ene rgy (V HE) gamma-rays 
was first pointed out by Nikishov ([1963). An original way 
of using ground-based TeV observations of distant sources 
to probe the CIB was given by Stecker et al. (1992) '). The 
fundamental idea is to look for absorption in the intrinsic 
spectrum as a result of electron-positron pair production 
by photon collisions 7TeV -|- 7cib — > e'^ + e~ . In such an 
interaction between a gamma-ray of energy (1 + z)E and 
an infrared photon of energy (1 + z)e, with E and e the 
observed energies at 2: = 0, the pair production thresh- 
old is Ee {I + zf {I - cose) > 2 (mc^)^ where 6 is the 



angle between photons and m the rest mass of t he ele c- 
tron. The cross-section can be written as (Heitler l|l960l )): 

a = C (1 - /3^) (2/3 (/3^ -2) + {3- f3^) In (j^)) cm^ with 

/3 = (l-2(mc2) /{Ee{l-cose)){l + zY\ and ( = 

1.25 X 10^^^. If the infrared photons have a density number 
n (e) de cm~^, the corresponding optical depth for attenua- 
tion is 



riE) 



Ho 



d (cos9) ■ 



de n{e)a{E,e,e) 



where et = 2(mc^)V(-E(l - cose){l + zf), z^ is 
the redshift of the source and c the speed of light. This 
expression is strictly valid for an Einstein-de Sitter universe 
but remains an accurate approximation in the ACDM 
paradigm for low redshifts. The detected flux is then atten- 
uated by a factor of e~'^^^\ The CIB energy distribution is 
assumed to be independent of z as the 7-ray source redshift 
is very low (0.034). The maximum cross-section is reached 
for an infrared photon wavelength of Xcib ~ Ac „ "^ y where 
Ac — h/ijnc) is the Compton wavelength of the electron. 
Therefore, 7-photons with energy between a few TeV and 
20 TeV "see" CIB photons with wavelengths between 3.5 
and 100 ^m. 

To allow for an intuitive understanding of the exponen- 
tial absorption (not only as a function of the quantity of the 
absorbing material, as it should be, but also as a function 
of the energy of the gamma-rays) , one can assume that the 
CIB is roughly constant, i.e. e^n{e) ~ cte. The optical depth 
for an energy kE can then be written as: 



T{kE) ex / d(cos6 



1 — cost 



feE(l+2)^{l 



(„c2) = 

~T7 



de 



a{kEe,9) 



By a change of variable e -^ ke, it is straightforward to see 
that T{kE) ^ k X t{E). As ^ = $0 x e"^^^', this explains 
why a high energy gamma-ray will be exponantially more 
absorbed than a low energy one. 



4 UNFOLDING METHOD AND PHYSICAL 
CRITERIA 

To turn the absorption corrected spectra into a lower limit 
on the Hubble parameter, the intrinsic spectral energy 
density of Mrk 501 is assumed to be concave {i.e. with a 
negative second derivative) in the multi-TeV region. This 
quite c ommo n approach has be en de veloped in Renault 
et al. (|200ll 'l and Guy et al. (|2000l '). The conservative 
hypothesis is based on the fact that no natural physical 
process can re-inject energy above the Inverse-Compton 
bump maximum. Both Klein-Nishina effect and auto- 
absorption within the source would only cause the flux to 
decrease more and more rapidly as a function of energy. 
This is obviously true for leptonic models (a good review 
of theoretical perspectives, where each model produces a 
concave spectrum, can be found in Sauge (2004)) but also 
for more complex hadronic models where this concave shape 
is also expected, either because of proton-initiated-cascade 
(often modeled as a broken power-law around 3 TeV), or 
due to the inclusion of /i-synchrotron radiation (Mannheim, 
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private communication). Even in the extreme case where 
proton synchrotron radi ation is at the origin of the TeV 
bump (Aharonian (2000)), the shape of the emission 
remains concave. Looking at the absorption corrected 
spectra, it cannot be totally excluded that the maximum of 
the so-called Inverse-Compton peak is not yet reached at 
20 TeV. Physical parameters (essentially magnetic field and 
Doppler factor) required to produce such a spectral energy 
distribution maximum above 17 TeV are substantially 
disfavored but, even in this case, the pF[v) shape should 
remain concave. This latter point is demonstra ted b y 
BeppoSAX 1-100 keV measurements (Plan et al. (|l998l ')') 
showing a clearly concave spectrum before the synchrotron 
bump maximum (around 100 keV) which is supposed to be 
mimicked by the TeV-spectrum before the so-called Inverse- 
Compton bump maximum. This hypothesis is reinforced 
by the fact that the sub- TeV slope, which is independent 
of the CIB density beyond 3.5 /im, effectively refiects the 
X-ray slope in the keV range. This behavior (together with 
the corre lated variability reported e.g. b y Aharonian et 
al. (|l99g ) and Djannati-Atai et al. (| 19991 )) indicates that 
the same population of particles is at the origin of both 
X-ray and 7-ray emissions, whatever this population is. 
In particular, the self-synchro-Compton model fits sa t- 
isfactorily the absorption corrected data (Guy et al. ([20001)). 
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Figure 2. Unfolded gamma-ray spectra and associated fits for 
(from bottom to top) : no CIB (equivalent to Hq — > oo), Hq = 
90 km.s-l.Mpc"^ and Hq = 50 km.s-^.Mpc"^. 



Furthermore, in addition to those "theoretically mo- 
tivated" arguments, observations show that the spectrum 
should indeed be concave. This can be checked by consid- 
ering AGNs with an IC bump seen at lower energies and 
therefore unaffected by the CIB absorption. This can easily 
be seen for 3C273 with the average spectrum compiled 
from 30 years of observations as reported in Tiirler et al. 
(|l99 9r). This is also true for Cen A (see, e.g., Chiaberge et 
al. ([2P01)) which is the nearest radiogalaxy (z=0.0018) and 
one of the best studied. Although with less significance this 
can even be checked f or M 87 (which is not a blazar), as 
shown in Lenain et al. (|2008l ). The FSRQ PK S 0521 -36 (see, 
e.g. the summary presented in Lenain et al. l|2008l ). mostly 
from Giommi et al. (|2002l ') and the NED), which is known 
for oscillating between a Seyfert-like and a BL Lac state, 
also exhibits a fully concave spectrum. This is also true for 
the blazar 3C 454. 3, as measured during the 2007 July flare 
(Ghisellini et al. ((2007])). During flares, good quality data 
reinforcing the concavity assumpti on ar e also available for 
PKS 2155-304 (Aharonian et al. (|2007l ')') and for 3 C 279 
(see, e.g., the summary plot in Bottcher et al. (|2007|)). 

It should be stressed that requiring the spectrum to 
be concave is "more demanding" than simply requiring 
the end of the TeV spectrum not to exhibit a too hard 
spectral index, as often used to derive limits on the CIB. 
Our approach evades the argument of Stec ker at al. \200Tl ) 
questioning the results of Aharonian et al. (|2006l ) . 

In order to quantify the concavity of the absorption 
corrected spectrum, a parabolic fit is performed from 6 
to 21 TeV in the plane {log{u),log{i'Fi')). This function, 
which is simple, is chosen for its constant second derivative 
a — d?{log{vFv))/d?(log{v)) which avoids the choice of 
a particular test-energy and fits satisfactorily the data. 
Over such a small interval it is obviously meaningful to 



assume the spectrum to be locally parabolic. The previous 
physical constraint on the concavity of the TeV spectrum 
simply reads as a < 0. As explained in the next section, 
the uncertainties are carefully taken into account. The 
parameter a is computed for different values of the Hubble 
parameter and the corresponding value is rejected if a is 
positive within the associated errors. Fig.[2]shows the result 
of the parabolic fit in several cases, superimposed with the 
experimental points. The fit under-estimates the energy 
density of the hardest photons, ensuring that the local 
second derivative is under-estimated, thus making the test 
conservative. 



5 RESULTS ON THE HUBBLE PARAMETER 

To derive statistically meaningful constraints on the Hubble 
parameter, errors on both the CIB and the measured TeV 
spectrum have been taken into account. As far as the 
gamma-ray spectrum is concerned, we have not only con- 
sidered the statistical uncertainties but also the systematic 
errors, as evaluated by the HEGRA team (Aharonian et al. 
(|l999bl )). The errors on the CIB density are not as easy to 
take into account as those on the AGN spectrum as their 
propagation is far from trivial. For this reason, we have 
performed a full Monte-Carlo simulation. Each of the eight 
considered CIB points between 0.3 and 160 /xm is randomly 
generated according to its central value and uncertainty, 
as displayed by the shaded region of Fig[T] Linear interpo- 
lations are then performed in the log-log plane to obtain 
a CIB realization. For each CIB realization, the unfolded 
spectrum is computed [i.e. the optical depth is numerically 
evaluated for the considered energies with the formula 
given in section 3) and fitted with a parabola between 6 
and 21 TeV. The error Aa on the second derivative a of 
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the parabola (due to uncertainties in the measurement of 
the gamma-ray spectrum) is computed according to three 
different statistical methods iterated up to the convergence 
point. The value of the Hubble parameter is rejected at 
xa if a — xAa > 0. The procedure is finally repeated 
to generate many (typically ~ 1000) realizations of the 
CIB. Confidence levels can be computed in this way by 
estimating the percentage of rejected realizations for each 
tested value of the Hubble parameter. It has been checked 
that our results are marginally dependent upon the number 
of points used in the fit and that the x^ P^r degree of 
freedom is, in each case, low enough to ensure a satisfactory 
fit quality. 

At the 68% confidence level, the lower limit de- 
rived on Ho, taking into account all the uncertainties is: 
Ho > 74 km.s~^.Mpc^^. This is the main result of this 
article. Although it should be underlined that even if a la 
limit remains rather weak on a statistical basis, this bound 
is meaningful -at this confidence level- when compared 
with other results. The history of the measured value of the 
Hubble parameter is very chaotic and different (incompati- 
ble) values have been taken as granted at different times. As 
reminded in t he rec ent review article "the Hubble constant" 
by Jackson ([ 20071 ), the current estimates range between 
60 and 75 km.s~'^.Mpc~^. In this context, our lower limit 
-ffo > 74 km.s~^.Mpc~^ is unquestionably relevant. 

As far as Cepheids are concerned, the main measure- 
ments are Ho = 73 km.s"\Mpc"^ (Riess et al. (2005)), 
with statistical errors of 4 km.s~^.Mpc~^ and system- 
atic errors of 5 km s~^ Mpc~^ on the o ne ha nd and 
Ho = 62.3 km.s~\Mpc"^ (Sandage et al. (|2006l ')'). with 
statistical errors of 1.3 km s~^ Mpc~^ and systematic errors 
of 5 km.s~^.Mpc~^ on the other hand. Our results clearly 
disfavor the Sandage estimate, even within its quoted 
uncertainty. As far as gravitational lenses are concerned, 
a sophisticat ed me ta-analysis has recently been performed 
(Oguri et al. (120071 )) using Monte-Carlo methods to account 
for quantities such as the presence of clusters around 
the main lens galaxy and the variation in profile slopes. 
The result obtained is Ho = 68 ± 6 ± 8 km.s"^Mpc"^ 
Meanwhile, an approach modelling simultaneously 10 of 
the 18 time-delay lenses with non-parametric models, has 
led (Saha et al. (2006 )) to Ho = 72 ± 12 km.s-\Mpc"\ 
Our results allow a significant la narrowing of those 
uncertainty intervals. As far as the Sunyaev-Zeldovich effect 
is con cerned, recent measurements (see, e.g. Jones et al. 
(|2005f) l are compatible with those estimates but the errors 
are too large to allow for any improvement. 

The most reliable measures available to date are 
probably those reported in the "final results from the 
Hubble Space Telescope Key Project" (Freedman et al. 
\2OOW 1). They concluded that the full analysis of all data 
gives Ho = 72 ± 8 km.s"^Mpc"\ Combining this HST 
allowed interval with our lower limit increases the mean 
value to Ho ~ 76 km.s~^.Mpc~^. However, the statistical 
significance of this estimate should be taken with care as 
the accurate distribution is not known. Nevertheless, our 
result Ho > 74 km.s~^.Mpc~^ clearly favors the upper end 
of the HST interval and substantially reduces the allowed 
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Figure 3. Lower limit at the 90% confidence level on tlic Hubble 
constant as a function of the rescaling factor of the CIB density. 
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than currently estimated. Each point is obtained by analysing 
1000 Monte-Carlo realizations of the CIB. 



parameter space. A less constraining, although interesting, 
result was obtained by Chandra (Bonamente et al. (l200y)) 
leading to Ho ^ 77 ± 12 km s"^ Mpc"^ with a larger 15 % 
uncertainty. Finally, WMA P has published in the five-years 
results (Komatsu et al. \200St i) an impressive estimate 
at 7^0 = 72 ± 2.6 km s~^ Mpc~^. It should, however, be 
strongly emphasized that WMAP results rely on many 
different priors and do not allow for an independent 
determination of Ho due to its degeneracy with the total 
curvature of the Universe. For example, every decrease of 
20 km.s~^.Mpc~^ in Ho increases the total density of the 
Universe by 0.1 in units of the closure density. The WMAP 
data by themselves, without any further assumptions or 
extra hypothesis, do not supply a significant constraint on 
Ho. 

If the statistical significance is to be improved to 90%, 
the lower limit derived in this work is significantly reduced to 
-ffo > 56 km.s~^.Mpc~^. This is, of course, due to large un- 
certainties on the CIB density. This broadening is somehow 
analogous to what happens for othe r app roaches to measure 
_ffo- For example, Tegmark et al. (2006) have shown that 
when relaxing the constraint of the equation of state of the 
dark energy, the 90% confidence level estimate using WMAP 
data is Ho G [61 - 84] km.s"\Mpc"^ 



CONCLUSION AND PROSPECTS 

The lower limit on _ffo derived in this letter is totally 
independent of the other available measurements. It is 
conservative, both due to the method which tends to 
under-estimate the high energy flux density and to the 
weak hypothesis regarding the spectrum concavity. Internal 
effects, such as Klein-Nishina cutoff and self-absorption, 
lead to very sharp multi-TeV intrinsic spectra. Furthermore, 
and most importantly, the CIB density is underestimated 
as only lower limits are used. Nevertheless, the errors 
on those lower limits have been exhaustively taken into 
account thanks to a Monte-Carlo simulation. Our result 
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allows a substantial la narrowing of previously obtained 
uncertainty intervals, favoring the higher possible values 
of Ho- It also closes the window on exotic "low Hubble 
constant" scenarios (Nugent at al. ( 1995J )) and on some 
alternatives to the cosmological concordance model (see, 
e.g., Blanchard et al. (|2003l ) where h ~ 0.46 is assumed.) 

Clearly, this work aims at providing a first hint in this 
direction and a first order of magnitude of what can be 
obtained by this method. Although the statistical signifi- 
cance of our results has been computed with care, several 
improvements could be expected in the future. First, other 
AGNs will pr obably be detected b y H.E .S.S (Aharonian 
et al. iioOMk^ MAGIC (Bastieri (|2006l )) or VERITAS 
(Krennrich ( 20071 ')) and sh ould b e included in the analysis. 
GLAST (Latronico et al. (|2007|)) data, at lower energies, 
could even be useful. This is very promising to account for 
more sources at different redshifts. However, the current 
situation for known AGNs is not expected to be improved 
significantly for a simple reason: the most important part 
of the spectrum used to constrain the Hubble parameter is 
the high-energy tail and larger telescopes will not improve 
the sensitivity in this range (just because the effective area 
is determined by the area of the shower and not by the 
area of the mirror). Then, improvements can be expected 
from measurements of the GIB. Fig |3] displays the evolution 
of the 2ct lower limit as a function of the GIB uncertainty 
scaling factor. A scaling factor x means that the errors are 
assumed to be x times smaller than currently estimated. 
The lower bound is, as expected, improved by smaller 
uncertainties and reaches a limit when the propagation 
of errors induced by the GIB uncertainties becomes much 
smaller than the effect of gamma-ray uncertainties associ- 
ated with the measurement of the AGN spectrum. For each 
point of the plot, 1000 Monte-Carlo realizations have been 
computed and analyzed according to the method described 
in the previous section. 

In the near future, measurements should indeed be im- 
proved between 60 and 11 /xm thanks to Herschel (see, e.g. 
Franceschini et al. l|2006r )) and, in the far future, between 
5 and 60 fj,m thanks to JWST (see, e.g. Windhorst et al. 
(1200a)). This could make this approach quite competitive. 

Finally, if gamma telescopes become so efficient that 
many more blazars are detected, each of them exhibiting 
a redshift-dependant cutoff in the spect rum, then the ap- 
proach suggested by Salamon et al. (119941 ) could be tractable 
and would lead not only to a bound but to a measurement 
of the Hubble parameter. However, to distinguish between 
a GIB absorption effect and an intrinsic cutoff, a very large 
number of data points would be required. 
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